Todo:

library(tidyverse)
library(readxl)
library(plotly)
ghg_emissions_clean <- ghg_emissions_data %>% 
  janitor::clean_names() %>% 
  mutate(units = "megatonnes of co2 equivelant") %>% 
  rename(emissions = emissions_mt_co2e) %>% 
  filter(emission_year != "BaseYear") %>% 
  mutate(emission_year = as.numeric(emission_year)) %>% 
  select(ccp_mapping, source_name, pollutant, year = emission_year, value = 
           emissions, units)
ghg_emissions_clean %>% 
  write_csv("data/clean_data/ghg_emissions.csv")
ghg_emissions_data %>% 
  names()
[1] "National Communication Categories" "SG Source Sector"                  "CCP mapping"                       "IPCC"                             
[5] "SourceName"                        "Pollutant"                         "EmissionYear"                      "Emissions (MtCO2e)"               
ghg_emissions_data %>% 
  select()
input <- list()

input$col_choice = "national_communication_categories"
ghg_emissions_clean %>%
  group_by_(input$col_choice, "emission_year") %>% 
  summarise(total_ghg_emissions = sum(emissions))
`summarise()` regrouping output by 'national_communication_categories' (override with `.groups` argument)
ghg_emissions_data %>% 
  distinct(`National Communication Categories`)
ghg_emissions_data %>% 
  filter(EmissionYear != "BaseYear") %>% 
  mutate(EmissionYear = as.numeric(EmissionYear)) %>% 
  group_by(EmissionYear) %>% 
  summarise(total_ghg_emissions = sum(`Emissions (MtCO2e)`)) %>% 
  ggplot() +
  aes(x = EmissionYear, y = total_ghg_emissions) +
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks = seq(1990,2020,5)) +
  ylim(0, 80) +
  theme(legend.position = 0) +
  theme_bw()
`summarise()` ungrouping output (override with `.groups` argument)

ghg_emissions_data %>% 
  distinct(`CCP mapping`)
ghg_emissions_data %>% 
  distinct(`National Communication Categories`)
ghg_emissions_data %>% 
  filter(EmissionYear != "BaseYear") %>% 
  mutate(EmissionYear = as.numeric(EmissionYear)) %>% 
  group_by(`CCP mapping`, EmissionYear) %>% 
  summarise(total_ghg_emissions = sum(`Emissions (MtCO2e)`)) %>% 
  ggplot() +
  aes(x = EmissionYear, y = total_ghg_emissions, group = `CCP mapping`, colour = `CCP mapping`) +
  geom_line() +
  geom_point() +
  scale_x_continuous(breaks = seq(1990,2020,5))
`summarise()` regrouping output by 'CCP mapping' (override with `.groups` argument)

ghg_emissions_data %>% 
  filter(EmissionYear != "BaseYear") %>% 
  mutate(EmissionYear = as.numeric(EmissionYear)) %>% 
  filter(`National Communication Categories` != `CCP mapping`) %>% 
  select(`National Communication Categories`, `CCP mapping`) %>% 
  unique()

category_id,category_name,subcategory_id,subcategory_name,year,emissions,emission

emissions_sankey <- emissions_data %>% 
  select(ccp_mapping, source_name, pollutant, emission_year, emissions, units)
filtered_df <- ghg_emissions_clean %>% 
  select(ccp_mapping, source_name, pollutant, emission_year, emissions, units) %>% 
  filter(pollutant == "CO2") %>% 
  filter(emission_year == "2005")
total_emissions_for_gas <- filtered_df %>% 
  summarise(sum(emissions)) %>% 
  pull()

total_emissions_by_category <- filtered_df %>% 
  select(-pollutant) %>% 
  group_by(ccp_mapping) %>% 
  summarise(cat_sum = sum(emissions), .groups = 'drop_last')
categories <- filtered_df %>% 
  distinct(ccp_mapping) %>% 
  pull()

n_categories <- length(categories)

sources <- filtered_df %>% 
  distinct(source_name) %>% 
  pull()

n_sources <- length(sources)
n_sources
[1] 159
node_names <- c("Total", categories, sources, "Other")

node_names_df <- data.frame("name" = node_names)

total_sankey_tibble <- total_emissions_by_category %>%
  mutate(total = "Total") %>% 
  mutate(total = match(total, node_names) -1) %>% 
  mutate(ccp_mapping = match(ccp_mapping, node_names) -1) %>% 
  select(source = total,
         target = ccp_mapping,
         value = cat_sum)

total_filtered_emissions <- total_sankey_tibble %>% 
  summarise(sum(value)) %>% 
  pull()

other_emissions <- total_emissions_for_gas - total_filtered_emissions

total_other_sankey_tibble <- tibble(
  "source" = c(0),
  "target" = (match("Other", node_names) -1),
  "value" = c(other_emissions)
)


sub_sankey_tibble <- filtered_df %>% 
  select(- c(units, pollutant, emission_year)) %>% 
  mutate(ccp_mapping = match(ccp_mapping, node_names) -1,
         source_name = match(source_name, node_names) -1)

names(sub_sankey_tibble) = c("source", "target", "value")

sankey_tibble <- total_sankey_tibble %>% 
  bind_rows(sub_sankey_tibble) %>% 
  bind_rows(total_other_sankey_tibble)

links_matrix <- data.frame(as.matrix(sankey_tibble, byrow = TRUE, ncols = 3))

# Add a 'group' column to each connection:
links <- links_matrix %>% 
  mutate(group = case_when(
    source == 0 ~ paste("type_", target, sep = ""),
    source!=0 ~ paste("type_", source, sep = "")
  ))

nodes <- node_names_df
# Add a 'group' column to each node.
# All of them in the same group to make them the same colour
nodes$group <- as.factor(c("my_unique_group"))

emissions <- list()

emissions$nodes <- nodes
emissions$links <- links
total_filtered_emissions <- total_sankey_tibble %>% 
  summarise(sum(value)) %>% 
  pull()

other_emissions <- total_emissions_for_gas - total_filtered_emissions

total_other_sankey_tibble <- tibble(
  "source" = c(0),
  "target" = (match("Other", node_names) -1),
  "value" = c(other_emissions)
)

sub_sankey_tibble <- filtered_tibble %>% 
  select(-category_id, -subcategory_id, -year) %>% 
  mutate(category_name = match(category_name, node_names) -1,
         subcategory_name = match(subcategory_name, node_names) -1)

names(sub_sankey_tibble) = c("source", "target", "value")

sankey_tibble <- total_sankey_tibble %>% 
  bind_rows(sub_sankey_tibble) %>% 
  bind_rows(total_other_sankey_tibble)

links_matrix <- data.frame(as.matrix(sankey_tibble, byrow = TRUE, ncols = 3))

# Add a 'group' column to each connection:
links <- links_matrix %>% 
  mutate(group = case_when(
    source == 0 ~ paste("type_", target, sep = ""),
    source!=0 ~ paste("type_", source, sep = "")
  ))

nodes <- node_names_df
# Add a 'group' column to each node.
# All of them in the same group to make them the same colour
nodes$group <- as.factor(c("my_unique_group"))

emissions <- list()

emissions$nodes <- nodes
emissions$links <- links
make_sankey_dfs <- function(data, userYear, userGas) {
  n_categories <- data %>% 
    distinct(category_name) %>% 
    nrow()
  
  total_emissions_for_gas <- data %>% 
    filter(emission == userGas()) %>% 
    filter(year == userYear()) %>%
    summarise(sum(emissions)) %>% 
    pull()
  
  filtered_tibble <- data %>%
    filter(emission == userGas()) %>% 
    select(-emission) %>% 
    filter(year == userYear()) %>% 
    filter(emissions > userResolution())
  
  total_emissions_by_cat <- filtered_tibble %>%
    group_by(category_name) %>% 
    summarise(cat_sum = sum(emissions), .groups = 'drop_last')
  
  categories <- filtered_tibble %>%
    distinct(category_name) %>% 
    pull()
  
  subcategories <- filtered_tibble %>%
    distinct(subcategory_name) %>% 
    pull()
  
  node_names <- c("Total", categories, subcategories, "Other")
  
  node_names_df <- data.frame("name" = node_names)
  
  total_sankey_tibble <- total_emissions_by_cat %>%
    mutate(total = "Total") %>% 
    mutate(total = match(total, node_names) -1) %>% 
    mutate(category_name = match(category_name, node_names) -1) %>% 
    select(source = total,
           target = category_name,
           value = cat_sum)
  
  total_filtered_emissions <- total_sankey_tibble %>% 
    summarise(sum(value)) %>% 
    pull()
  
  other_emissions <- total_emissions_for_gas - total_filtered_emissions
  
  total_other_sankey_tibble <- tibble(
    "source" = c(0),
    "target" = (match("Other", node_names) -1),
    "value" = c(other_emissions)
  )
  
  sub_sankey_tibble <- filtered_tibble %>% 
    select(-category_id, -subcategory_id, -year) %>% 
    mutate(category_name = match(category_name, node_names) -1,
           subcategory_name = match(subcategory_name, node_names) -1)
  
  names(sub_sankey_tibble) = c("source", "target", "value")
  
  sankey_tibble <- total_sankey_tibble %>% 
    bind_rows(sub_sankey_tibble) %>% 
    bind_rows(total_other_sankey_tibble)
  
  links_matrix <- data.frame(as.matrix(sankey_tibble, byrow = TRUE, ncols = 3))
  
  # Add a 'group' column to each connection:
  links <- links_matrix %>% 
    mutate(group = case_when(
      source == 0 ~ paste("type_", target, sep = ""),
      source!=0 ~ paste("type_", source, sep = "")
    ))
  
  nodes <- node_names_df
  # Add a 'group' column to each node.
  # All of them in the same group to make them the same colour
  nodes$group <- as.factor(c("my_unique_group"))
  
  emissions <- list()
  
  emissions$nodes <- nodes
  emissions$links <- links
  
  return(emissions)
}
make_sankey_dfs(emissions_sankey, userYear = 2005, userGas = "CH4", userResolution = 50)
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KIyBUb2RvOg0KDQotIGFkZCBjb2x1bW5zIGZvciB0b3RhbHMgZm9yIGVhY2ggcG9sbHV0YW50L3NlY3Rvcg0KDQpgYGB7cn0NCmxpYnJhcnkodGlkeXZlcnNlKQ0KbGlicmFyeShyZWFkeGwpDQpsaWJyYXJ5KHBsb3RseSkNCmBgYA0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfZGF0YSA8LSByZWFkX3hsc3goImRhdGEvcmF3X2RhdGEvc2NvdHRpc2gtZ2hnLWRhdGFzZXQtMjAxOC54bHN4IiwgMikNCmBgYA0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfY2xlYW4gPC0gZ2hnX2VtaXNzaW9uc19kYXRhICU+JSANCiAgamFuaXRvcjo6Y2xlYW5fbmFtZXMoKSAlPiUgDQogIG11dGF0ZSh1bml0cyA9ICJtZWdhdG9ubmVzIG9mIGNvMiBlcXVpdmVsYW50IikgJT4lIA0KICByZW5hbWUoZW1pc3Npb25zID0gZW1pc3Npb25zX210X2NvMmUpICU+JSANCiAgZmlsdGVyKGVtaXNzaW9uX3llYXIgIT0gIkJhc2VZZWFyIikgJT4lIA0KICBtdXRhdGUoZW1pc3Npb25feWVhciA9IGFzLm51bWVyaWMoZW1pc3Npb25feWVhcikpICU+JSANCiAgc2VsZWN0KGNjcF9tYXBwaW5nLCBzb3VyY2VfbmFtZSwgcG9sbHV0YW50LCB5ZWFyID0gZW1pc3Npb25feWVhciwgdmFsdWUgPSANCiAgICAgICAgICAgZW1pc3Npb25zLCB1bml0cykNCmBgYA0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfY2xlYW4gJT4lIA0KICB3cml0ZV9jc3YoImRhdGEvY2xlYW5fZGF0YS9naGdfZW1pc3Npb25zLmNzdiIpDQpgYGANCg0KDQpgYGB7cn0NCnBsb3RseTo6Z2dwbG90bHkoDQogIGdoZ19lbWlzc2lvbnNfY2xlYW4gJT4lIA0KICAgIGZpbHRlcih5ZWFyID09IDIwMTgpICU+JQ0KICAgIGZpbHRlcihwb2xsdXRhbnQgJWluJSBjKCJDTzIiLCAiTjJPIiwgIkNINCIpKSAlPiUgDQogICAgZ2dwbG90KCkgKw0KICAgIGFlcyh4ID0gZmFjdG9yKGNjcF9tYXBwaW5nLCBsZXZlbHMgPSByZXYobGV2ZWxzKGZhY3RvcihjY3BfbWFwcGluZykpKSksDQogICAgICAgIHkgPSB2YWx1ZSwgZmlsbCA9IHNvdXJjZV9uYW1lLA0KICAgICAgICB0ZXh0ID0gcGFzdGUwKCc8L2JyPiBTZWN0b3I6ICcsIGNjcF9tYXBwaW5nLA0KICAgICAgICAgICAgICAgICAgICAgICc8L2JyPiBFbWlzc2lvbnM6ICcsIHZhbHVlLA0KICAgICAgICAgICAgICAgICAgICAgICc8L2JyPiBTb3VyY2UgTmFtZTogJywgc291cmNlX25hbWUpKSArDQogICAgZ2VvbV9jb2wocG9zaXRpb24gPSAic3RhY2siKSArDQogICAgdGhlbWVfYncoKSArDQogICAgZmFjZXRfd3JhcCh+cG9sbHV0YW50KSArDQogICAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm5vbmUiKSArDQogICAgbGFicyh4ID0gIlNlY3RvciIsDQogICAgICAgICB5ID0gcGFzdGUwKCJFbWlzc2lvbnMgKCIsIGdoZ19lbWlzc2lvbnNfY2xlYW4kdW5pdHNbMV0sICIpIikpICsNCiAgICBjb29yZF9mbGlwKCksDQogICAgdG9vbHRpcCA9ICd0ZXh0Jw0KICApDQoNCmBgYA0KDQoNCmBgYHtyfQ0KZ2hnX2VtaXNzaW9uc19kYXRhICU+JSANCiAgbmFtZXMoKQ0KYGBgDQoNCmBgYHtyfQ0KZ2hnX2VtaXNzaW9uc19kYXRhICU+JSANCiAgc2VsZWN0KCkNCmBgYA0KDQoNCg0KYGBge3J9DQppbnB1dCA8LSBsaXN0KCkNCg0KaW5wdXQkY29sX2Nob2ljZSA9ICJuYXRpb25hbF9jb21tdW5pY2F0aW9uX2NhdGVnb3JpZXMiDQpgYGANCg0KDQoNCmBgYHtyfQ0KZ2hnX2VtaXNzaW9uc19jbGVhbiAlPiUNCiAgZ3JvdXBfYnlfKGlucHV0JGNvbF9jaG9pY2UsICJlbWlzc2lvbl95ZWFyIikgJT4lIA0KICBzdW1tYXJpc2UodG90YWxfZ2hnX2VtaXNzaW9ucyA9IHN1bShlbWlzc2lvbnMpKSAlPiUgDQogIGdncGxvdCgpICsNCiAgYWVzKHggPSBFbWlzc2lvblllYXIsIHkgPSB0b3RhbF9naGdfZW1pc3Npb25zLCBncm91cCA9IGBOYXRpb25hbCBDb21tdW5pY2F0aW9uIENhdGVnb3JpZXNgLCBjb2xvdXIgPSBgTmF0aW9uYWwgQ29tbXVuaWNhdGlvbiBDYXRlZ29yaWVzYCkgKw0KICBnZW9tX2xpbmUoKSArDQogIGdlb21fcG9pbnQoKSArDQogIHNjYWxlX3hfY29udGludW91cyhicmVha3MgPSBzZXEoMTk5MCwyMDIwLDUpKSArDQogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9IDApDQpgYGANCg0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGEgJT4lIA0KICBkaXN0aW5jdChgTmF0aW9uYWwgQ29tbXVuaWNhdGlvbiBDYXRlZ29yaWVzYCkNCmBgYA0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGEgJT4lIA0KICBmaWx0ZXIoRW1pc3Npb25ZZWFyICE9ICJCYXNlWWVhciIpICU+JSANCiAgbXV0YXRlKEVtaXNzaW9uWWVhciA9IGFzLm51bWVyaWMoRW1pc3Npb25ZZWFyKSkgJT4lIA0KICBncm91cF9ieShFbWlzc2lvblllYXIpICU+JSANCiAgc3VtbWFyaXNlKHRvdGFsX2doZ19lbWlzc2lvbnMgPSBzdW0oYEVtaXNzaW9ucyAoTXRDTzJlKWApKSAlPiUgDQogIGdncGxvdCgpICsNCiAgYWVzKHggPSBFbWlzc2lvblllYXIsIHkgPSB0b3RhbF9naGdfZW1pc3Npb25zKSArDQogIGdlb21fbGluZSgpICsNCiAgZ2VvbV9wb2ludCgpICsNCiAgc2NhbGVfeF9jb250aW51b3VzKGJyZWFrcyA9IHNlcSgxOTkwLDIwMjAsNSkpICsNCiAgeWxpbSgwLCA4MCkgKw0KICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAwKSArDQogIHRoZW1lX2J3KCkNCmBgYA0KDQoNCmBgYHtyfQ0KZ2hnX2VtaXNzaW9uc19kYXRhICU+JSANCiAgZGlzdGluY3QoYENDUCBtYXBwaW5nYCkNCmBgYA0KDQpgYGB7cn0NCg0KYGBgDQoNCg0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGEgJT4lIA0KICBkaXN0aW5jdChgTmF0aW9uYWwgQ29tbXVuaWNhdGlvbiBDYXRlZ29yaWVzYCkNCmBgYA0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGENCmBgYA0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2RhdGEgJT4lIA0KICBmaWx0ZXIoRW1pc3Npb25ZZWFyICE9ICJCYXNlWWVhciIpICU+JSANCiAgbXV0YXRlKEVtaXNzaW9uWWVhciA9IGFzLm51bWVyaWMoRW1pc3Npb25ZZWFyKSkgJT4lIA0KICBncm91cF9ieShgQ0NQIG1hcHBpbmdgLCBFbWlzc2lvblllYXIpICU+JSANCiAgc3VtbWFyaXNlKHRvdGFsX2doZ19lbWlzc2lvbnMgPSBzdW0oYEVtaXNzaW9ucyAoTXRDTzJlKWApKSAlPiUgDQogIGdncGxvdCgpICsNCiAgYWVzKHggPSBFbWlzc2lvblllYXIsIHkgPSB0b3RhbF9naGdfZW1pc3Npb25zLCBncm91cCA9IGBDQ1AgbWFwcGluZ2AsIGNvbG91ciA9IGBDQ1AgbWFwcGluZ2ApICsNCiAgZ2VvbV9saW5lKCkgKw0KICBnZW9tX3BvaW50KCkgKw0KICBzY2FsZV94X2NvbnRpbnVvdXMoYnJlYWtzID0gc2VxKDE5OTAsMjAyMCw1KSkNCmBgYA0KDQpgYGB7cn0NCmdoZ19lbWlzc2lvbnNfZGF0YSAlPiUgDQogIGZpbHRlcihFbWlzc2lvblllYXIgIT0gIkJhc2VZZWFyIikgJT4lIA0KICBtdXRhdGUoRW1pc3Npb25ZZWFyID0gYXMubnVtZXJpYyhFbWlzc2lvblllYXIpKSAlPiUgDQogIGZpbHRlcihgTmF0aW9uYWwgQ29tbXVuaWNhdGlvbiBDYXRlZ29yaWVzYCAhPSBgQ0NQIG1hcHBpbmdgKSAlPiUgDQogIHNlbGVjdChgTmF0aW9uYWwgQ29tbXVuaWNhdGlvbiBDYXRlZ29yaWVzYCwgYENDUCBtYXBwaW5nYCkgJT4lIA0KICB1bmlxdWUoKQ0KYGBgDQpjYXRlZ29yeV9pZCxjYXRlZ29yeV9uYW1lLHN1YmNhdGVnb3J5X2lkLHN1YmNhdGVnb3J5X25hbWUseWVhcixlbWlzc2lvbnMsZW1pc3Npb24NCg0KYGBge3J9DQplbWlzc2lvbnNfc2Fua2V5IDwtIGVtaXNzaW9uc19kYXRhICU+JSANCiAgc2VsZWN0KGNjcF9tYXBwaW5nLCBzb3VyY2VfbmFtZSwgcG9sbHV0YW50LCBlbWlzc2lvbl95ZWFyLCBlbWlzc2lvbnMsIHVuaXRzKQ0KYGBgDQoNCg0KYGBge3J9DQpnaGdfZW1pc3Npb25zX2NsZWFuDQpgYGANCg0KDQpgYGB7cn0NCmZpbHRlcmVkX2RmIDwtIGdoZ19lbWlzc2lvbnNfY2xlYW4gJT4lIA0KICBzZWxlY3QoY2NwX21hcHBpbmcsIHNvdXJjZV9uYW1lLCBwb2xsdXRhbnQsIGVtaXNzaW9uX3llYXIsIGVtaXNzaW9ucywgdW5pdHMpICU+JSANCiAgZmlsdGVyKHBvbGx1dGFudCA9PSAiQ08yIikgJT4lIA0KICBmaWx0ZXIoZW1pc3Npb25feWVhciA9PSAiMjAwNSIpDQpgYGANCg0KYGBge3J9DQp0b3RhbF9lbWlzc2lvbnNfZm9yX2dhcyA8LSBmaWx0ZXJlZF9kZiAlPiUgDQogIHN1bW1hcmlzZShzdW0oZW1pc3Npb25zKSkgJT4lIA0KICBwdWxsKCkNCg0KdG90YWxfZW1pc3Npb25zX2J5X2NhdGVnb3J5IDwtIGZpbHRlcmVkX2RmICU+JSANCiAgc2VsZWN0KC1wb2xsdXRhbnQpICU+JSANCiAgZ3JvdXBfYnkoY2NwX21hcHBpbmcpICU+JSANCiAgc3VtbWFyaXNlKGNhdF9zdW0gPSBzdW0oZW1pc3Npb25zKSwgLmdyb3VwcyA9ICdkcm9wX2xhc3QnKQ0KYGBgDQoNCmBgYHtyfQ0KY2F0ZWdvcmllcyA8LSBmaWx0ZXJlZF9kZiAlPiUgDQogIGRpc3RpbmN0KGNjcF9tYXBwaW5nKSAlPiUgDQogIHB1bGwoKQ0KDQpuX2NhdGVnb3JpZXMgPC0gbGVuZ3RoKGNhdGVnb3JpZXMpDQoNCnNvdXJjZXMgPC0gZmlsdGVyZWRfZGYgJT4lIA0KICBkaXN0aW5jdChzb3VyY2VfbmFtZSkgJT4lIA0KICBwdWxsKCkNCg0Kbl9zb3VyY2VzIDwtIGxlbmd0aChzb3VyY2VzKQ0KYGBgDQoNCmBgYHtyfQ0Kbl9zb3VyY2VzDQpgYGANCg0KDQpgYGB7cn0NCm5vZGVfbmFtZXMgPC0gYygiVG90YWwiLCBjYXRlZ29yaWVzLCBzb3VyY2VzLCAiT3RoZXIiKQ0KDQpub2RlX25hbWVzX2RmIDwtIGRhdGEuZnJhbWUoIm5hbWUiID0gbm9kZV9uYW1lcykNCg0KdG90YWxfc2Fua2V5X3RpYmJsZSA8LSB0b3RhbF9lbWlzc2lvbnNfYnlfY2F0ZWdvcnkgJT4lDQogIG11dGF0ZSh0b3RhbCA9ICJUb3RhbCIpICU+JSANCiAgbXV0YXRlKHRvdGFsID0gbWF0Y2godG90YWwsIG5vZGVfbmFtZXMpIC0xKSAlPiUgDQogIG11dGF0ZShjY3BfbWFwcGluZyA9IG1hdGNoKGNjcF9tYXBwaW5nLCBub2RlX25hbWVzKSAtMSkgJT4lIA0KICBzZWxlY3Qoc291cmNlID0gdG90YWwsDQogICAgICAgICB0YXJnZXQgPSBjY3BfbWFwcGluZywNCiAgICAgICAgIHZhbHVlID0gY2F0X3N1bSkNCg0KdG90YWxfZmlsdGVyZWRfZW1pc3Npb25zIDwtIHRvdGFsX3NhbmtleV90aWJibGUgJT4lIA0KICBzdW1tYXJpc2Uoc3VtKHZhbHVlKSkgJT4lIA0KICBwdWxsKCkNCg0Kb3RoZXJfZW1pc3Npb25zIDwtIHRvdGFsX2VtaXNzaW9uc19mb3JfZ2FzIC0gdG90YWxfZmlsdGVyZWRfZW1pc3Npb25zDQoNCnRvdGFsX290aGVyX3NhbmtleV90aWJibGUgPC0gdGliYmxlKA0KICAic291cmNlIiA9IGMoMCksDQogICJ0YXJnZXQiID0gKG1hdGNoKCJPdGhlciIsIG5vZGVfbmFtZXMpIC0xKSwNCiAgInZhbHVlIiA9IGMob3RoZXJfZW1pc3Npb25zKQ0KKQ0KDQoNCnN1Yl9zYW5rZXlfdGliYmxlIDwtIGZpbHRlcmVkX2RmICU+JSANCiAgc2VsZWN0KC0gYyh1bml0cywgcG9sbHV0YW50LCBlbWlzc2lvbl95ZWFyKSkgJT4lIA0KICBtdXRhdGUoY2NwX21hcHBpbmcgPSBtYXRjaChjY3BfbWFwcGluZywgbm9kZV9uYW1lcykgLTEsDQogICAgICAgICBzb3VyY2VfbmFtZSA9IG1hdGNoKHNvdXJjZV9uYW1lLCBub2RlX25hbWVzKSAtMSkNCg0KbmFtZXMoc3ViX3NhbmtleV90aWJibGUpID0gYygic291cmNlIiwgInRhcmdldCIsICJ2YWx1ZSIpDQoNCnNhbmtleV90aWJibGUgPC0gdG90YWxfc2Fua2V5X3RpYmJsZSAlPiUgDQogIGJpbmRfcm93cyhzdWJfc2Fua2V5X3RpYmJsZSkgJT4lIA0KICBiaW5kX3Jvd3ModG90YWxfb3RoZXJfc2Fua2V5X3RpYmJsZSkNCg0KbGlua3NfbWF0cml4IDwtIGRhdGEuZnJhbWUoYXMubWF0cml4KHNhbmtleV90aWJibGUsIGJ5cm93ID0gVFJVRSwgbmNvbHMgPSAzKSkNCg0KIyBBZGQgYSAnZ3JvdXAnIGNvbHVtbiB0byBlYWNoIGNvbm5lY3Rpb246DQpsaW5rcyA8LSBsaW5rc19tYXRyaXggJT4lIA0KICBtdXRhdGUoZ3JvdXAgPSBjYXNlX3doZW4oDQogICAgc291cmNlID09IDAgfiBwYXN0ZSgidHlwZV8iLCB0YXJnZXQsIHNlcCA9ICIiKSwNCiAgICBzb3VyY2UhPTAgfiBwYXN0ZSgidHlwZV8iLCBzb3VyY2UsIHNlcCA9ICIiKQ0KICApKQ0KDQpub2RlcyA8LSBub2RlX25hbWVzX2RmDQojIEFkZCBhICdncm91cCcgY29sdW1uIHRvIGVhY2ggbm9kZS4NCiMgQWxsIG9mIHRoZW0gaW4gdGhlIHNhbWUgZ3JvdXAgdG8gbWFrZSB0aGVtIHRoZSBzYW1lIGNvbG91cg0Kbm9kZXMkZ3JvdXAgPC0gYXMuZmFjdG9yKGMoIm15X3VuaXF1ZV9ncm91cCIpKQ0KDQplbWlzc2lvbnMgPC0gbGlzdCgpDQoNCmVtaXNzaW9ucyRub2RlcyA8LSBub2Rlcw0KZW1pc3Npb25zJGxpbmtzIDwtIGxpbmtzDQoNCg0KYGBgDQoNCg0KDQoNCg0KDQpgYGB7cn0NCnRvdGFsX2ZpbHRlcmVkX2VtaXNzaW9ucyA8LSB0b3RhbF9zYW5rZXlfdGliYmxlICU+JSANCiAgc3VtbWFyaXNlKHN1bSh2YWx1ZSkpICU+JSANCiAgcHVsbCgpDQoNCm90aGVyX2VtaXNzaW9ucyA8LSB0b3RhbF9lbWlzc2lvbnNfZm9yX2dhcyAtIHRvdGFsX2ZpbHRlcmVkX2VtaXNzaW9ucw0KDQp0b3RhbF9vdGhlcl9zYW5rZXlfdGliYmxlIDwtIHRpYmJsZSgNCiAgInNvdXJjZSIgPSBjKDApLA0KICAidGFyZ2V0IiA9IChtYXRjaCgiT3RoZXIiLCBub2RlX25hbWVzKSAtMSksDQogICJ2YWx1ZSIgPSBjKG90aGVyX2VtaXNzaW9ucykNCikNCg0Kc3ViX3NhbmtleV90aWJibGUgPC0gZmlsdGVyZWRfdGliYmxlICU+JSANCiAgc2VsZWN0KC1jYXRlZ29yeV9pZCwgLXN1YmNhdGVnb3J5X2lkLCAteWVhcikgJT4lIA0KICBtdXRhdGUoY2F0ZWdvcnlfbmFtZSA9IG1hdGNoKGNhdGVnb3J5X25hbWUsIG5vZGVfbmFtZXMpIC0xLA0KICAgICAgICAgc3ViY2F0ZWdvcnlfbmFtZSA9IG1hdGNoKHN1YmNhdGVnb3J5X25hbWUsIG5vZGVfbmFtZXMpIC0xKQ0KDQpuYW1lcyhzdWJfc2Fua2V5X3RpYmJsZSkgPSBjKCJzb3VyY2UiLCAidGFyZ2V0IiwgInZhbHVlIikNCg0Kc2Fua2V5X3RpYmJsZSA8LSB0b3RhbF9zYW5rZXlfdGliYmxlICU+JSANCiAgYmluZF9yb3dzKHN1Yl9zYW5rZXlfdGliYmxlKSAlPiUgDQogIGJpbmRfcm93cyh0b3RhbF9vdGhlcl9zYW5rZXlfdGliYmxlKQ0KDQpsaW5rc19tYXRyaXggPC0gZGF0YS5mcmFtZShhcy5tYXRyaXgoc2Fua2V5X3RpYmJsZSwgYnlyb3cgPSBUUlVFLCBuY29scyA9IDMpKQ0KDQojIEFkZCBhICdncm91cCcgY29sdW1uIHRvIGVhY2ggY29ubmVjdGlvbjoNCmxpbmtzIDwtIGxpbmtzX21hdHJpeCAlPiUgDQogIG11dGF0ZShncm91cCA9IGNhc2Vfd2hlbigNCiAgICBzb3VyY2UgPT0gMCB+IHBhc3RlKCJ0eXBlXyIsIHRhcmdldCwgc2VwID0gIiIpLA0KICAgIHNvdXJjZSE9MCB+IHBhc3RlKCJ0eXBlXyIsIHNvdXJjZSwgc2VwID0gIiIpDQogICkpDQoNCm5vZGVzIDwtIG5vZGVfbmFtZXNfZGYNCiMgQWRkIGEgJ2dyb3VwJyBjb2x1bW4gdG8gZWFjaCBub2RlLg0KIyBBbGwgb2YgdGhlbSBpbiB0aGUgc2FtZSBncm91cCB0byBtYWtlIHRoZW0gdGhlIHNhbWUgY29sb3VyDQpub2RlcyRncm91cCA8LSBhcy5mYWN0b3IoYygibXlfdW5pcXVlX2dyb3VwIikpDQoNCmVtaXNzaW9ucyA8LSBsaXN0KCkNCg0KZW1pc3Npb25zJG5vZGVzIDwtIG5vZGVzDQplbWlzc2lvbnMkbGlua3MgPC0gbGlua3MNCmBgYA0KDQoNCg0KDQoNCg0KDQoNCmBgYHtyfQ0KbWFrZV9zYW5rZXlfZGZzIDwtIGZ1bmN0aW9uKGRhdGEsIHVzZXJZZWFyLCB1c2VyR2FzKSB7DQogIG5fY2F0ZWdvcmllcyA8LSBkYXRhICU+JSANCiAgICBkaXN0aW5jdChjYXRlZ29yeV9uYW1lKSAlPiUgDQogICAgbnJvdygpDQogIA0KICB0b3RhbF9lbWlzc2lvbnNfZm9yX2dhcyA8LSBkYXRhICU+JSANCiAgICBmaWx0ZXIoZW1pc3Npb24gPT0gdXNlckdhcygpKSAlPiUgDQogICAgZmlsdGVyKHllYXIgPT0gdXNlclllYXIoKSkgJT4lDQogICAgc3VtbWFyaXNlKHN1bShlbWlzc2lvbnMpKSAlPiUgDQogICAgcHVsbCgpDQogIA0KICBmaWx0ZXJlZF90aWJibGUgPC0gZGF0YSAlPiUNCiAgICBmaWx0ZXIoZW1pc3Npb24gPT0gdXNlckdhcygpKSAlPiUgDQogICAgc2VsZWN0KC1lbWlzc2lvbikgJT4lIA0KICAgIGZpbHRlcih5ZWFyID09IHVzZXJZZWFyKCkpICU+JSANCiAgICBmaWx0ZXIoZW1pc3Npb25zID4gdXNlclJlc29sdXRpb24oKSkNCiAgDQogIHRvdGFsX2VtaXNzaW9uc19ieV9jYXQgPC0gZmlsdGVyZWRfdGliYmxlICU+JQ0KICAgIGdyb3VwX2J5KGNhdGVnb3J5X25hbWUpICU+JSANCiAgICBzdW1tYXJpc2UoY2F0X3N1bSA9IHN1bShlbWlzc2lvbnMpLCAuZ3JvdXBzID0gJ2Ryb3BfbGFzdCcpDQogIA0KICBjYXRlZ29yaWVzIDwtIGZpbHRlcmVkX3RpYmJsZSAlPiUNCiAgICBkaXN0aW5jdChjYXRlZ29yeV9uYW1lKSAlPiUgDQogICAgcHVsbCgpDQogIA0KICBzdWJjYXRlZ29yaWVzIDwtIGZpbHRlcmVkX3RpYmJsZSAlPiUNCiAgICBkaXN0aW5jdChzdWJjYXRlZ29yeV9uYW1lKSAlPiUgDQogICAgcHVsbCgpDQogIA0KICBub2RlX25hbWVzIDwtIGMoIlRvdGFsIiwgY2F0ZWdvcmllcywgc3ViY2F0ZWdvcmllcywgIk90aGVyIikNCiAgDQogIG5vZGVfbmFtZXNfZGYgPC0gZGF0YS5mcmFtZSgibmFtZSIgPSBub2RlX25hbWVzKQ0KICANCiAgdG90YWxfc2Fua2V5X3RpYmJsZSA8LSB0b3RhbF9lbWlzc2lvbnNfYnlfY2F0ICU+JQ0KICAgIG11dGF0ZSh0b3RhbCA9ICJUb3RhbCIpICU+JSANCiAgICBtdXRhdGUodG90YWwgPSBtYXRjaCh0b3RhbCwgbm9kZV9uYW1lcykgLTEpICU+JSANCiAgICBtdXRhdGUoY2F0ZWdvcnlfbmFtZSA9IG1hdGNoKGNhdGVnb3J5X25hbWUsIG5vZGVfbmFtZXMpIC0xKSAlPiUgDQogICAgc2VsZWN0KHNvdXJjZSA9IHRvdGFsLA0KICAgICAgICAgICB0YXJnZXQgPSBjYXRlZ29yeV9uYW1lLA0KICAgICAgICAgICB2YWx1ZSA9IGNhdF9zdW0pDQogIA0KICB0b3RhbF9maWx0ZXJlZF9lbWlzc2lvbnMgPC0gdG90YWxfc2Fua2V5X3RpYmJsZSAlPiUgDQogICAgc3VtbWFyaXNlKHN1bSh2YWx1ZSkpICU+JSANCiAgICBwdWxsKCkNCiAgDQogIG90aGVyX2VtaXNzaW9ucyA8LSB0b3RhbF9lbWlzc2lvbnNfZm9yX2dhcyAtIHRvdGFsX2ZpbHRlcmVkX2VtaXNzaW9ucw0KICANCiAgdG90YWxfb3RoZXJfc2Fua2V5X3RpYmJsZSA8LSB0aWJibGUoDQogICAgInNvdXJjZSIgPSBjKDApLA0KICAgICJ0YXJnZXQiID0gKG1hdGNoKCJPdGhlciIsIG5vZGVfbmFtZXMpIC0xKSwNCiAgICAidmFsdWUiID0gYyhvdGhlcl9lbWlzc2lvbnMpDQogICkNCiAgDQogIHN1Yl9zYW5rZXlfdGliYmxlIDwtIGZpbHRlcmVkX3RpYmJsZSAlPiUgDQogICAgc2VsZWN0KC1jYXRlZ29yeV9pZCwgLXN1YmNhdGVnb3J5X2lkLCAteWVhcikgJT4lIA0KICAgIG11dGF0ZShjYXRlZ29yeV9uYW1lID0gbWF0Y2goY2F0ZWdvcnlfbmFtZSwgbm9kZV9uYW1lcykgLTEsDQogICAgICAgICAgIHN1YmNhdGVnb3J5X25hbWUgPSBtYXRjaChzdWJjYXRlZ29yeV9uYW1lLCBub2RlX25hbWVzKSAtMSkNCiAgDQogIG5hbWVzKHN1Yl9zYW5rZXlfdGliYmxlKSA9IGMoInNvdXJjZSIsICJ0YXJnZXQiLCAidmFsdWUiKQ0KICANCiAgc2Fua2V5X3RpYmJsZSA8LSB0b3RhbF9zYW5rZXlfdGliYmxlICU+JSANCiAgICBiaW5kX3Jvd3Moc3ViX3NhbmtleV90aWJibGUpICU+JSANCiAgICBiaW5kX3Jvd3ModG90YWxfb3RoZXJfc2Fua2V5X3RpYmJsZSkNCiAgDQogIGxpbmtzX21hdHJpeCA8LSBkYXRhLmZyYW1lKGFzLm1hdHJpeChzYW5rZXlfdGliYmxlLCBieXJvdyA9IFRSVUUsIG5jb2xzID0gMykpDQogIA0KICAjIEFkZCBhICdncm91cCcgY29sdW1uIHRvIGVhY2ggY29ubmVjdGlvbjoNCiAgbGlua3MgPC0gbGlua3NfbWF0cml4ICU+JSANCiAgICBtdXRhdGUoZ3JvdXAgPSBjYXNlX3doZW4oDQogICAgICBzb3VyY2UgPT0gMCB+IHBhc3RlKCJ0eXBlXyIsIHRhcmdldCwgc2VwID0gIiIpLA0KICAgICAgc291cmNlIT0wIH4gcGFzdGUoInR5cGVfIiwgc291cmNlLCBzZXAgPSAiIikNCiAgICApKQ0KICANCiAgbm9kZXMgPC0gbm9kZV9uYW1lc19kZg0KICAjIEFkZCBhICdncm91cCcgY29sdW1uIHRvIGVhY2ggbm9kZS4NCiAgIyBBbGwgb2YgdGhlbSBpbiB0aGUgc2FtZSBncm91cCB0byBtYWtlIHRoZW0gdGhlIHNhbWUgY29sb3VyDQogIG5vZGVzJGdyb3VwIDwtIGFzLmZhY3RvcihjKCJteV91bmlxdWVfZ3JvdXAiKSkNCiAgDQogIGVtaXNzaW9ucyA8LSBsaXN0KCkNCiAgDQogIGVtaXNzaW9ucyRub2RlcyA8LSBub2Rlcw0KICBlbWlzc2lvbnMkbGlua3MgPC0gbGlua3MNCiAgDQogIHJldHVybihlbWlzc2lvbnMpDQp9DQpgYGANCg0KDQpgYGB7cn0NCm1ha2Vfc2Fua2V5X2RmcyhlbWlzc2lvbnNfc2Fua2V5LCB1c2VyWWVhciA9IDIwMDUsIHVzZXJHYXMgPSAiQ0g0IiwgdXNlclJlc29sdXRpb24gPSA1MCkNCmBgYA0KYGBge3J9DQoNCmBgYA0KDQo=